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The thermodynamics of supersymmetric Yang-Mills theories is studied by computing the 
two-loop correction to the canonical free energy and to the equation of state for theories 
with 16, 8 and 4 supercharges in any dimension 4 < d < 10, and in two dimensions at finite 
volume. In the four-dimensional case we also evaluate the first non-analytic contribution 
in the 't Hooft coupling to the free energy, arising from the resummation of ring diagrams. 
To conclude, we discuss some applications to the study of the Hagedorn transition in string 
theory in the context of Matrix strings and speculate on the possible physical meaning of 
the transition. 
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1. Introduction 



The study of the thermal properties of supersymmetric Yang-Mills (SYM) and super- 
string theories has received a boost @ @ H @ @ after the D-brane revolution @ . More 
recently, the conjecture by Maldacena of a duality between string theory on anti-de Sitter 
(AdS) backgrounds and the large-N limit of SYM theories || has further motivated the 
study of these issues in the hope that it could lead to a clarification of the mechanisms of 
confinement in non-abelian gauge theories ||10|| . 

Beyond its purely cosmological/phenomeno logical relevance, the study of SYM ther- 
modynamics finds interesting applications in the study of near-extremal black holes ||11|| 



and D-branes ||12]| [13]. When looking at the effective theory of coincident D-branes there 
are, generically, two regimes associated with the values of the effective gauge couplings 
9sQ (Q being the number of coincident D-branes or the charges of the black hole). When 
g s Q < 1 the field theory limit is well described by a perturbative SYM theory living on 
the world-volume of the D-brane. On the other hand, taking g s Q > 1 the weakly coupled 
D-brane picture is not appropriate any more. Using Maldacena's conjecture, however, it is 
possible to relate the field theory limit in this non-perturbative regime (large 't Hooft cou- 
pling) with a supergravity computation on some background of the form AdSd x (Spheres). 
Corrections to the leading strong coupling result of order 0(l/g s Q) are then associated 
with higher dimensional terms (a/ corrections) in the supergravity effective action. 

The study of the corrections to the leading results on both sides (weak |14| and strong 
|l5|][|l(| coupling) has been done for the conformal M = 4 SYM theory in d = 4 and a 
tendency in both curves to meet was detected (see however [|I7|]). In this note we will 
try to achieve a twofold objective. First, to compute the two-loop free energy of SYM 
theories with 16, 8 and 4 supercharges in various dimensions. The first class of theories 
are specially interesting because of their potential application to non-conformal versions of 
the Maldacena conjecture @. Second, to obtain the first non-analytic correction in the 't 
Hooft coupling, of order 0[(gYMN) 3 / 2 ], to the two-loop result in four dimensional theories 
arising from the resummation of ring diagrams. We will analyze also with some detail 
the two-dimensional case, where the infrared divergences will be handled by putting the 
system at finite volume. 

The physics of the high temperature string gas has been a recurrent issue in string 
theory (for a sample of papers from the "golden age" of thermal strings see (18|] [|19|] ) . In 
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recent years we have gained interesting insights about the physical meaning of the Hagedorn 
divergence, in spite of the fact that a full and detailed understanding of the problem seems 
to be still at large. Although this will not be the main subject of this note, we will try to 
discuss some aspects of the Hagedorn transition that could be enlightened by our results 
on SYM thermodynamics. In particular we will use our study of the two-dimensional 
SYM theories to try to get some qualitative information about the Hagedorn transition 
using Matrix strings as a non-perturbative definition of Type-IIA superstrings. In spite of 
being non-conclusive, we hope the discussion to be helpful in shedding some light to such 
a confusing issue. 

The present paper is organized as follows: in the next Section, the two-loop corrections 
to the thermal free energy of SYM theories with 16, 8 and 4 supercharges will be computed 
in any dimension 4 < d < 10 using dimensional reduction from the corresponding maximal 
M = 1 SYM theory. In Section 2.3 we compute the next correction to the two-loop free 
energy for SYM theories in four dimensions. Section 2.4 will be devoted to the study 
of the two-dimensional case at finite volume. Finally, in Section 3 we will summarize 
the conclusions and discuss some possible application of our results to the study of the 
Hagedorn transition in Matrix string theory. 



2. Two-loop free energy of SYM theories 

2.1. Super •symmetric Yang-Mills theories in various dimensions 

In this section we will compute the next-to-leading contribution to the canonical free 
energy of supersymmetric Yang-Mills theories with 16, 8 and 4 supercharges. In order to 
keep the analysis general, we will start with M = 1 SYM in T> dimensions, whose dynamics 



is governed by the action [£0 

S = ! d°x Tr 



4 #M 



F AB F AB +i4,T A D A ^ 



where A, B = 0, . . . , T> — 1 and both the gauge fields and spinors are in the adjoint repre- 
sentation of U(N). We get theories with different number of supercharges by choosing the 
appropriate value of T> for which that number #sc is maximal: 



#sc~- 


= 16 - 


^ ^max 


= 10 


#sc 


= 8- 




= 6 


#sc 


= 4- 


^max 


= 4. 
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In addition, different conditions on fermions have to be imposed in order to keep the number 
of bosonic and fermionic degrees of freedom equal. Thus, when P max = 10 fermions have to 
be taken Majorana-Weyl, while for D max = 6 and D max = 4 they satisfy Weyl conditions 
(actually in d = 4 we can choose the fermions to be either Majorana or Weyl, both 
conditions being equivalent |2G[|). This ensures that the number of physical bosonic and 
fermionic degrees of freedom will be equal to P max — 2. 

In general, however, we will be interested in SYM theories with ifsc supercharges 
in dimensions d < £> max . This theories can be obtained by dimensional reduction of the 
corresponding maximal Af = 1 SYM theory in D = P max • Thus, we can parameterize 
any (i-dimensional SYM theory with any number of supercharges by specifying both d and 
the maximal dimension £> max from which it is obtained by dimensional reduction. In this 
way, starting with £> max = 10 (Af = 1 in d = 10) we get Af = 1 in d = 8, Af = 2 in d = 6, 
Af = 4 in d = 4 and Af = 8 in d = 2. Starting instead with P max = 6 (Af = 1 in d = 6) we 
will have Af = 2 in d = 4 and Af = 4 in d = 2. Finally, if we take D max = 4 we can retrieve 
A" = 2 in d = 2. (for odd dimensions 2n — 1 we have the Af corresponding to dimension 
2n). 



2.2. Two-loop free energy for d > 4 

With this in mind, we can proceed to compute the canonical free energy in pertur- 
bation theory for any supersymmetric Yang-Mills theory characterized by (T> max ,d), by 
writing down the contribution of vacuum Feynman diagrams of Af = 1 in T> = £> max SYM 
and restricting internal momentum in loops to d dimensions. That way we are able to 
keep track of the contribution of gauge bosons and scalars (as well as their supersymmet- 
ric partners) without having to consider a larger number of diagramsi. The final result, 
of course, will depend on (d, P ma x)- 

As a warmup exercise, we will compute the one-loop free energy density. In the 
maximal Af = 1 theory we have three diagrams, a bosonic loop, a fermionic loop and the 



2 For example, using this trick one can get the result of ref. [14] by computing, instead of ten, 
only four two-loop Feynman diagrams. 



ghost loop, which after multiplying by their corresponding degeneracy factor respectively 
give (we use the notation of ref. |2TJ) 

H3)i-ioop = JiV 2 D ma jf [d d P]\ogP 2 - JiV 2 (X? max -2)rf [d d P]logP 2 
z T(p) ^ 7{-P} 

— N 2 j£ [d d P] log P 2 
T(p) 

where by (P) and {P} we represent, respectively, bosonic (periodic) and fermionic (anti- 
periodic) boundary conditions along the Euclidean time and the factor of iV 2 is due to 
the fact that all fields are in the adjoint representation of U(N). After a straightforward 
computation we arrive at 

[d d P] logP 2 ^ W 7 f^ r log(p 2 + u 2 n ) = A - ^PcW d 
T(p) P , , J (27rr 7T 2 

jf [d d P] logP 2 = 1 £ / 7 f^ r log(p 2 + c 2 ) = A + - 2^)^)/^. 

T{P} P (27T) 01 7T2 

Ao is a regularized vacuum energy that will cancel after summing all contributions. The 
total result for the one-loop free energy is thus 

H3)i-ioo P = -^^C(«Q(1 - 2- d )(2V ax - 2)N 2 (3- d - (2.1) 

7T2 

Next we get the two-loop corrections to this result. Thus, we must sum the contribu- 
tions of the four Feynman diagrams of Fig. 1. corresponding to Af = 1 SYM in T> = D max 
where internal momenta is restricted to a <i-dimensional space-time (one of whose direc- 
tions is the compactified Euclidean time). Proceeding this way and after some algebra we 
find the contribution of each independent diagram (using the Feynman-'t Hooft gauge) 

n = j^ew -2)ivi J - 2 {i P) lddp] -k) {i P) lddp] -k 
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Fig. 1: Feynman diagrams contributing to the two- loop canonical free en- 
ergy. Solid lines represent fermions, wavy lines gauge bosons and dashed lines 
Fadeev-Popov ghosts. 

Here Tr 1 is the dimension of the spinors of the maximal SYM theory which in all the 
cases under study (Af = 1 in "D max = 10, 6, 4) equals £> ma x — 2. Using this fact we can add 
all the above contributions and find the following expression for the two-loop free energy 
density 

W)2-ioo P = \g 2 Y M N\V m ^ - 2) 2 (jf {p ld d P}^2 - jf {p} ld d P]^J ■ (2-2) 

The integrals appearing between brackets contain both the zero and finite temperature 
part of J 7 (/3)2-ioop- They can be easily computed by performing the sum first, to give 



(P) [a J P2 J (2*)*-i2u, p + 



4 



d d ~ x p 


N p 


(27r) d " 1 


u p 




n p 


(2 7 r) d " 1 


UJ p 



(2.3) 



<{P} 

where u p = p and N p , n p are the Bose-Einstein and Fermi-Dirac distribution functions 
respectively 

N p = e ^ p _ l > n P = e /3u p + l ■ 



The first thing to be said about eq. (|2.2|) is that the zero temperature (ultraviolet 
divergent) contribution cancels out between the bosonic and the fermionic integral. This 
is just a consequence of supersymmetry since the vacuum energy should not be corrected 
at zero temperature if supersymmetry is to be preserved by the vacuum |J22| . This is the 
reason why this cancellation occurs not only for the conformal Af = 4, SYM4 (P ma x = 10) 



Tfl , but for all SYM theories under study. 

Notice that although all SYM theories in dimension higher than 4 are non- 
renormalizable, the two-loop finite temperature free energy is well behaved in the ultravi- 
olet. This is a consequence of the fact that the ultraviolet region in the thermal integrals 
is effectively cut-off for momenta p ^> T and therefore the temperature dependent part of 
the amplitudes is, to a great extent, insensitive to ultraviolet ambiguities. Of course, di- 
vergences should have been taken care of in the zero temperature sector by an appropriate 
cut-off in momenta A (although some protected observables, like the vacuum energy, will 
be finite due to supersymmetry). In that case, consistency will require T < A. 

To get an analytical expression for the free energy we can evaluate the integrals ap- 
pearing in ( p.3|) for generic values of the dimension! 



(27T) d - 1 LO p r(^fi) ^ ' ^ ' 

[ J^PIhL = (1 _ 2 3 - d ) 22 ~y~Z C(d - 2)T(d - 2)p 2 ~ d 7 
so after substituting in (|2.3|) and (|2.2|) we get 



2-d 



^ r (/5)2-ioo P = 9ymN" 



2 

(v max - 2) 22 ;_! C J1~ ? (i - 2*- d )r(d - 2 



7T 2 



4-2d 



Actually, the fermionic integral can be written in terms of the bosonic one using 

2 1 1 



which is just a realization of the well-known relation between the one-loop free energy of a bosonic 
and a fermionic quantum field, Ff er (/3) = i 7 b os (/3) — 2Fb os (2/3) 
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/3 4 " d + OM M Nf 



which together with the one-loop contribution can be written as 
F{(3) = -Y 2 /T d /?EW2) c(d)(1 _ 2 - d )(V max - 2) 

2d n 2 

(2.4) 

Incidentally, the two- loop correction to the free energy is always positive for d > 1, so it 
always tends to increase the (negative) one-loop contribution. 

From ( |2.4| ) we see explicitly how corrections to the one-loop result come in powers of 
the 't Hooft coupling gyM-N which is kept fixed in the large-N limit. Since in d-dimensions 
the Yang-Mills coupling constant gyM nas dimension of (Energy) 4- ^ the condition for 
the perturbative expansion to be reliable is the effective dimensionless coupling at a given 
temperature to be small 

g 2 eff = gy M Nf3 4 - d <^l. 

Notice that for d > 4, for which the SYM theory is non-renormalizable, perturbative 
corrections to the free energy will be governed by a small parameter in the low temperature 
limit, corresponding to the fact that the theory is well behaved in the infrared. Actually, 
higher order corrections to formula ( |2.4| ) have better and better infrared behavior as we 
increase the order in perturbation theory. At the same time the ultraviolet divergences 
worsen in the zero temperature sector, while in the temperature-dependent part of the 
amplitude the ultraviolet behavior is smoother due to the presence of Boltzmann factors 
that effectively cut-off momenta beyond a scale of orderS T. 

From the two-loop canonical free energy we can obtain the corrections to the equation 
of state of SYMrf. We first compute the canonical entropy density S(T) as a function of 
the temperature and invert it up to terms of order 0[(gyMN) 2 ] to get 



d-3 

S N 



4 The better ultraviolet behavior of the temperature-dependent sector of the theory does not 
guarantee its finiteness in higher loops; for example, the free energy in SYMio is ultraviolet 
divergent at three loops as can be seen by thermal averaging the one-loop effective action, which 
contains a F 4 term that scales quadratically with the ultraviolet cut-off. The resulting thermal 
averaged divergent part is of order 0[(g Y M Y) 2 ] [fL6|l , as corresponds to a three-loop contribution. 
I thank A. Tseytlin for pointing this out to me. 
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where the numerical coefficients JF , T\ are defined from the free energy ( |2.4|) by T = 
-T N 2 T d + {g YM N)N 2 T 1 T 2d - /L . Now we can substitute into the internal energy density 
£ = T + TS with the result 



2.3. Thed = 4 case 



d 
d-1 



2 \ ( ^ 

1 + i9YMN) Jd^Wo\N^d 



d-4 
d-1 



+ 0[{g 2 YM Nf] 



We will now concentrate our attention in the four-dimensional case. Taking d = 4 in 
we find 

.2 



•F(/3)d=4 = -^iV 2 /r 4 (D max - 2) 

o 



y-^Pmax-2)^ M 7V 



(2.5) 



As a check, we can evaluate (|2.4|) for the superconformal A/" = 4 SYM4 theory, which is 
obtained by dimensional reduction of J\f = 1 in SYM10 (i.e. V max = 10) to give 



•F(/3)at=4 = -A 2 /3 



2 «-4 



7T 



6 



1 



which indeed agrees with the result of ref. [[14| . In the equation of state, the loop correction 
just results in a renormalization of the overall numerical factor (this also happens in the 
non-conformal cases V max = 6, 4) 



n- 



2tt s 



(1 + 



Generically, the next contribution to ( p.4|) is naively given by three-loop diagrams of 
order 0[(gY M N) 2 }. However in four dimensions, as it happens in QCD [p4]|||21|| ) at three 
loop level there are already uncanceled infrared divergences that have to be cured by sum- 
ming over ring diagrams. This gives a non-analytic (of order 0[(<7y M iV) 3 / 2 ]) contribution 
to the free energy, representing a mild failure of perturbation theory due to the infrared 
ambiguities. The evaluation of this term is essentially equivalent to dressing the Aq and 
scalar propagators in loops by introducing the effect of Debye screening and thermal mass 
for the scalars. To leading order in the 't Hooft coupling, the electric (Debye) mass can be 
easily computed from the static limit of the one-loop self-energy to give 

1 



m^ = limffi a (0,p) = -(P r 
p— »o 4 



2)g' YM NT* + 0[(g YM NY 



(2.6) 
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while for the scalars we have 

1 



,2 



limn fla (0,p) = ~(V max -2)g YM NT z + 0[(g YM N) 2 ]. (2.7) 

p^O o 

In order to get the 0[(g YM N) 3 / 2 ] terms in the free energy we use the technique of 



ref. | [2 If , and rewrite the original Lagrangian density as 

£sym 4 = (csYM^ + lmliTrAlSp^o + ^mlf^^] - i m^Tr^X.o - \™\ ^Tr0f 

\ i=l / i=l 

(2.8) 

where (pi are the n s = T> max — 4 adjoint scalars in the theory and the electric mass only 
affects to the zero-frequency component of the Aq field (cf. |[24|| ). The strategy now is to 
treat the last two terms as a perturbation to the Lagrangian density between brackets. This 
results in a reorganization of perturbation theory in which the ring-diagram contribution 
can be easily evaluated. 

The first thing will be to compute again the one-loop free energy density, including 
now the effect of the masses in the Lagrangian ( p.8|) and, at the same time, adding new 
one-loop diagrams containing vertices associated with the counterterms. Expanding the 
results up to order 0[(g Y m-^) 2 ] we find! 

HP)\-r op = -W)i-ioo P + ^N 2 T [ml + (P max - 4)m 3 J + 0[(g YM N) 2 } (2.9) 

with J? r (/9)i-ioop given by (|2.1| ). Proceeding similarly with the two-loops diagrams of Fig. 
1, we get 

^(/?)™p = HPh-ioo P ~ ^N 2 T [ml + (X? max - A)m%] + 0[(g YM N) 2 }. 

So we are left with the following final result for the "2^ -loop" contribution to the free 
energy density 

W) 2 ± -ioo P = -^ N " T Hi + - 4)to 3 J (2.10) 

where the values of the thermal masses are given by eqs. ( |2.6|) and (|2.7| ). It is important 
to notice that this term is always negative for all 4 < P max < 10. 



5 The trick of dimensional reduction is no longer useful here because the thermal mass dis- 
tinguishes between scalar and gauge boson propagators. Thus we have to compute all diagrams 
separately. 
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The only thing left now will be to add ( 2.10 ) to the two-loop result (|2.5| ). In particular, 
doing so for the superconformal M = 4, SYM4 and evaluating the numerical coefficients, 
we find (A 2 = g 2 YM N) 



T{(3) N=A = -iV 2 T 4 [1.645 - 0.250A 2 + 0.234A 3 + £>(A 4 )]. 

Next terms in the perturbative expansion in four dimensions will be of order O (A 4 ) and 
C(A 4 logA) and can be also evaluated using the strategy employed in [^T[ (or up to order 



0(X 5 ) using |^5[ |2~5|] ). However, finite temperature perturbation theory is expected to break 
down at order C[(<7ym-^0 3 ] ]27j. Since this failure of perturbation theory is associated with 
the infrared sector of the theory, supersymmetry is not expected to solve the problem or 
even improve the situation. As in QCD ||28||||26|| , some kind of non-perturbative analysis 
will be needed in order to compute higher orders. Actually, the general structure of the 
series in A is important in trying to decide whether there is a phase transition occurring 
at some intermediate value of the 't Hooft coupling that precludes the extrapolation of 
supergravity physics into the gauge theory domain ]Tj 



2-4- SYM2 thermodynamics on S\ x R 

When d < 3 the analysis of the quantum corrections to the one loop free energy gets ad- 
ditional complications due to the hard infrared divergences that afflict super-renormalizable 
theories. For d = 3, we see that expression ( |2.4|) diverges because of a £(1) factor. In prin- 
ciple this can be cured, as usual, by computing the thermal masses and inserting them into 
the propagators, thus regularizing the low-momentum behavior of the Feynman integrals. 
However, in the three-dimensional case the computation of the electric mass has to be done 
with extra care, since the one-loop corrections to the propagators are already infrared di- 
vergent. Thus, the electric mass has to be evaluated self-consistently a la Hartree-Fock 
||29|| . Anyway, we will not dwell in this case any further. 

The two-dimensional case, on the other hand, is more interesting from several points 
of view. The one that will concern us here is that M = 8, SYM2 describes the world- 



volume dynamics of Matrix strings |3Cj , a non-perturbative definition of Type-IIA super- 
strings. Naively, ( |2.4j ) is ill-defined for d = 2 due to the endemic infrared divergences of 
low-dimensional field theories. There are several ways in which this divergence can be 
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regularized. Here we will get rid of the problem by putting the system in finite bo xi of 
length L = 2nR. We will assume that the thermal wavelengths of the fundamental fields 
are much smaller than the global length of the box (3 <C L and restrict our analysis to the 
sector without Wilson lines (the "long strings" that characterize the matrix string phase). 
Once this is done, the only change in the computation of Feynman diagrams is that con- 
tinuous space momentum is discretized in units of 1/R and the momentum integrals have 
to be replaced by discrete sums 



(2.11) 



where the second sum is, as usual, over integer or half-integer Matsubara frequencies 
depending on the bosonic or fermionic character of the propagating field. 

In the one-loop approximation the relevant bosonic and fermionic determinants have 



been already computed in |Hj and the resulting one-loop free energy density can be cast 
in terms of modular functions 



T{&, L)i-ioo P = -^N 2 (V max - 2) lof 



9 a 



~ --(v r 



2)N 2 (3~ 



In computing the bosonic determinant, and in order to keep the argument of the logarithm 
dimensionless, we have added a /3-independent counterterm. In the infinite volume limit 
L — > oo we recover the one-loop result obtained in Sec. 2.2. 

Let us now go to the two-loop case. To compute the contribution to the free en- 
ergy density we can use formula (|2.2| ) provided we substitute the integration by the sum 



according to (|2.11| ). After making so, we find 



A-K 2 m 2 4n 2 n 2 
^ + 



i 



(2.12) 




Jj 



where we have made use of the Epstein zeta-function [52| and $(m, n) represents the 
quadratic form between the square brackets. It is interesting to notice here that the 



6 Actually, if we take the d — > 2 limit in expression ( ^.4|) we get a finite result with a two- 
loop correction independent of the temperature. However, since dimensional regularization is not 
reliable in dealing with infrared divergences we will not follow this procedure. 
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resulting regular zeta-function arises as the difference of two singular zeta-functions with 
£(l)-like divergences which cancel out. This is again due to the non-renormalization of the 
vacuum energy for super symmetric theories. 

Actually, the zeta-function in (|2.12j ) can be written itself in terms of ordinary modular 
functions (see first article in ||32|| ), so at the end we can write 



— loops 



1 



4tt 2 
N 2 



g$ M N 3 (V T 



2) 2 log 2 
L 2 



ft 



;2.i3) 



in such a way that one arrives at the following expression for the (l+2)-loop free energy 
density 



F([3,L) = - 1 ^N 2 \ (P max -2)log 



04 



l (3 



178 (if 



Atx- 



2Y(g^ M N)(L(3) log^ 



9? 



l 2/3 



~ "T^max 
/3<L 4 



2)N 2 P~ 



1 



1447T 



+ 0[{g 2 YM Nf 



2)g 2 YM NL' 



According to this formula, the natural effective dimensionless coupling in the large L limit is 
now 5^ = (gy M N)L 2 . The analysis with be reliable when 1 > (gy M N)L 2 > (gy M N)(3 2 . 
Again, the equation of state in the two dimensional case can be computed when /3 C L, 
with the result 

max 



£ 



nN 2 (V max -2) 



1 + 



1447T 



-QymNL 2 



(2.14) 



3. Conclusions and outlook: Hagedorn transition from SYM thermodynamics? 

In the present paper, the thermodynamics of supersymmetric Yang-Mills theories with 
16, 8 and 4 supercharges was studied in any dimension d > 4. We computed the two-loop 
correction to the free energy for these theories and found that it always has opposite sign 
to the leading (negative) one-loop result. In the four-dimensional case we also evaluated 
the correction to the free energy arising from the resummation of the ring diagrams, using 
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the technique of ref. [21], and found it to be negative. For lower dimensional (d < 3) SYM 
theories, the computation is plagued with infrared divergences that have to be regularized 
somehow. We studied in detail the two-dimensional case at finite volume (to regularize 
these infrared divergences) in the high temperature limit. Again we found a positive 
two-loop correction which scales as T 2 with an effective dimensionless coupling given by 
9 2 y M NL 2 . 

Before closing, let us make some remarks on the potential use of SYM thermodynam- 
ics in clarifying the issue of the Hagedorn transition. On general grounds, one can expect 
two possible resolutions to the Hagedorn problem: either non-perturbative effects drive 
the critical temperature to a maximum reachable temperature for the system or new fun- 
damental degrees of freedom appear at high energies, thus providing a picture for a phase 
transition (or a smooth crossover, depending on the details of the dynamics). Although at 
present there are no clear evidences as to which one of the two alternatives is physically 
realized in string/M-theory, some results HH and our still incomplete knowledge of the 
theory seem to hint in the direction of the second one. 

D-instanton corrections to the ther mo dynamical potentials have been studied in || 
with the result that they do not modify the critical behavior at the Hagedorn temperature. 
More recently, the authors of refs. || have included non-perturbative semiclassical ingredi- 
ents in the analysis of the physics of the Hagedorn transition at finite volume through the 
Horowitz-Polchinski correspondence principle |33[ , getting a picture in which the Hagedorn 
phase is bounded at high energies by a black hole phase. A similar situation occurs for a 
string gas on AdS backgrounds where, in the canonical ensemble, the Hagedorn transition 
is "screened" by the formation of an AdS black hole ||. 

A second approach to the problem would start with a non-perturbative formulation 
of string theory in terms of M-theoretic degrees of freedom, as it has been proposed in 
|f. Let us momentarily adhere ourselves to this latter path and, starting with the non- 
perturbative definition of the Type-IIA superstring provided by Matrix strings JJ(J , study 
the world-volume thermodynamics of Type-IIA strings in the micro canonical ensemble. 
The world-volume theory is governed by M = 8 SYM2 with the Yang-Mills coupling 
constant given by g\ M = l/(g 2 a'), with g s the string coupling constant. On the other 
hand, free field configurations are determined by the overall scale a . In the infrared, 
E <C gYM, the physics is dominated by "long string" excitations along the flat directions. 
It is in this regime in which Matrix strings reproduce, in the large-N limit, the multi-string 
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Type-IIA ensemble [00]. If the energy is increased, the system will begin to be excited 
along non-flat directions as well. At energies E 3> Qym the potential terms in the M = 8 
SYM2 theory will behave as a small perturbation and the system will enter a perturbative 
regime. Thermodynamics there is well defined, as we have seen from the previous analysis. 

It is tempting to try to make some connection between these two world-sheet regimes 
and the low/high energy regimes in the target space string theory. At low energies we 
have perturbative Type-IIA super string theory that, in ten open space-time dimensions, 
we know is characterized at high energies by a negative specific heat phase. This nega- 
tive specific heat phase is viewed as a breakdown of equipartition in energy, in the sense 
that most of the energy of the string ensemble is stored into one (or a small number) of 
highly excited strings |18[][]T9[]. From the philosophy of M-theory it seems quite reason- 
able to expect that if too much energy is stored into a single string some transition to 
non-perturbative (maybe eleven-dimensional) physics should take place, putting an end to 
the negative specific phase. Alternatively, a black hole could be formed before the system 
leaves the string regime [f|. In any case, the final conclusion would be that the Hagedorn 
phase will be bounded by a new phase into which the system will decay either via a smooth 
crossover or a phase transition. 

In the case at hand, however, it is not clear how to connect the world-volume theory 
with some kind of target picture. One of the difficulties lies in the fact that Matrix strings 
are formulated in the light-cone gauge, in which the space-time interpretation is rather 
obscure. Nonetheless, one can naively argue that the negative specific heat phase at inter- 
mediate energies o/ _1//2 < E <C g~ l a'~ 1 / 2 is bounded at high energy E ^> (7~ 1 a'~ 1 / 2 by a 
new phase with regular thermodynamics (i.e. positive specific heat) effectively described 
by a perturbative two-dimensional U(N) supersymetric Yang-Mills theory with sixteen 
supercharges in the large-N limit. If this were so, the transition between the low energy 
string phase and the new high energy phase would be through a first order phase transition 
across the unstable (negative specific heat) phase (cf. Carlitz in [T§]). The critical points 
would be determined by the Maxwell rule for the entropy, provided the complete profile of 
the micro canonical temperature T(E) is known. 

The space-time interpretation of such a phase is far from being straightforward. In 
the SYM2 perturbative regime (or directly in the free limit gy M ~ (g^a') -1 — > 0) the 
two-dimensional action is that of sigma-model in a "non-commutative" target space with 
matrix coordinates G Adj[U(N)]. Whether this indicates that the Hagedorn transition 
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corresponds!! to the nucleation of non-commutative bubbles in a commutative space-time 
is something that it is difficult to decide with our present knowledge of the theory. One of 
the problems to be clarified will be, for example, how the target space volume dependence 
of the extensive quantities emerges as a function of N. In any case, we stress that this 
extrapolation of the world-sheet picture to space-time physics is very speculative, and 
should be tested by a detailed computation. We hope to report on this elsewhere. 

In a sense, this picture can be regarded as dual to the one proposed in [|J. There, 
the Hagedorn transition is linked to the condensation of DO-branes and their low-energy 
dynamics will be U(N) super quantum mechanics with sixteen supercharges. Both descrip- 
tions could in principle be related by performing a T-duality along the ninth dimension 
and interchanging its role with the M-theory circle. 
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Note added 

After this paper appeared in the LANL hep-th archive, I learned directly from S.-J. 
Rey of his parallel and independent work with C. Kim on SYM thermodynamics, part of 
which overlaps with the results presented here and that has later appeared in [|53| . I would 
like also to thank A. Nieto and A. Tseytlin for their interesting remarks on the first version 
of the article. 



At least in those cases in which it is not preceded by the formation of black holes due to the 
corresponding principle. Actually, we can tune the string coupling constant g s , the volume and 
the total energy in such a way that the system avoid the correspondence line and thus we prevent 
the formation of black holes. 
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